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Ph! Abstract 

We study the distribution of lepton pairs from the second hghtest 
, neutrahno decay X2 ~* followed by / — > y^l. The distribution of 

the ratio of lepton transverse momenta At shows peak structure if 
mu < is required. The peak position A^^"^ is described by 

a simple function of the ino and slepton masses in the mu ~ limit. 
5^ \ When a moderate mu cut is applied, A^^^ depends on the X2 velocity 

distribution, but the dependence would be corrected by studying the 



X 



lepton Pt distribution. A^^^ and the edge of mu distributions are 
used to determine the mass parameters involved in the decay for pa- 
rameters of interest to LHC experiments. For some cases the ino and 
slepton masses may be determined within 10% by the lepton distribu- 
tion only independent of model assumptions. Correct combinations of 
A^^^ and m'lf^'^ would be identified even if different X2 decay chains 
are co-existing. The analysis could be extended to the Tevatron energy 
scale or other cascade decays. 



1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) [||] is one of the most 
promising extensions of the Standard Model. It offers a natural solution of 
the hierarchy problem, amazing gauge coupling unification, and dark matter 
candidates. If Nature chooses low energy supersymmetry (SUSY), sparticles 
will be found for sure, as they will be copiously produced at future colliders 
such as Large Hadron Collider (LHC) at CERN or TeV scale e~^e~ linear 
colliders (LC) proposed by DESY, KEK, and SLAG. LHC would be a great 
discovery machine. Squarks and gluinos with mass less than a few TeV would 
be found unless the decay patterns are non-canonical 0. 

On the other hand, the MSSM suffers severe flavor changing neutral cur- 
rent (FCNC) constraints if no mass relation is imposed on sfermion mass 
parameters [Q. Various proposals have been made for the mechanism to in- 
corporate SUSY breaking in "our sector" , trying to offer natural explanations 
of such mass relations @]. In short, it would be very surprising if sparticles 
are found in future collider experiments — The discovery is not the final 
goal, but it is the beginning of a new quest for "the mechanism" of SUSY 
breaking. 

Measurements of soft breaking masses would be an important aspect of 
the study of SUSY, because different SUSY breaking mechanisms predict 
different sparticle mass patterns. Studies at the Tevatron and LHC would 
suffer from substantial uncertainties and backgrounds compared to an LC, 
such as luminosity error, combinatorial backgrounds, and unknown initial 
energy. While the discovery of sparticles is guaranteed at the LHC, detailed 
studies there would be challenging. Therefore it is very interesting to see the 
ultimate precision of supersymmetric studies at the LHC. 
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It is possible to determine masses of sparticles from the measurement of 
end points of invariant mass distributions [Q, ^ |^ . For the minimal super- 
gravity (MSUGRA) and gauge mediated (GM) models, there was substantial 
success for the parameter points where the decay of the second lightest neu- 
tralino to lepton pair X2 ~^ ^^Xi is detected with substantial statistics. For 
some case, one would be able to not only determine all MSUGRA parame- 
ters, but also to measure the masses of some sparticles, using the edges and 
end points of invariant mass distributions involving jets and leptons. The 
systematic errors of such analyses may be controlled if the acceptance near 
the end points and (jet) energy resolution are known. 

Detailed studies in this direction have been performed, and we do not 
repeat these here. In this paper, we instead study the ratio of lepton Pj- 
(lepton Pt asymmetry At = Pt2/ ^tv Pt2 < Pti ) the decay X2 ~^ 
llXi- The information has been used in previous analyses ^, ^ in the context 
of global fits of MSUGRA parameters. We show that it is possible to make 
a direct connection between the peak structure of the asymmetry A^"^^ and 
the ratio of the lepton energies in the neutralino rest frame Ae hj using 
events with mu < We also point out that systematics due to the 

X2 velocity distribution would be small and reduced further if one includes 
the Pt distribution of the hardest lepton in the fit. Using the mu end point 
and the peak position of the At distribution, one can at least determine two 
degrees of freedom of the three parameters involved in the X2 decay, m^o, 
m^o and mj. The measurements are based on lepton distributions only and 
free from uncertainty due to jet energy smearing. 

The organization of this paper is as follows. In section 2, we analyze the 
MSUGRA points which were studied in 0, |[, where squark and gluino decays 
are the dominant sources of X2- We concentrate on the case where X2 ~^ is 
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open and followed by I — > Ixi- We find that the At distribution has a peak 
if mil ~ is required. In the limit where niu ~ 0, the peak necessarily 

agrees with the ratio of lepton energies = E12/ En in the X2 ^^st frame for 
any value of the X2 velocity. is a simple function of the ino and slepton 
masses. We show that a small mu cut promises smaller systematic errors 
by comparing distributions for different neutralino velocities. In section 3, 
we show Monte Carlo simulations for the MSUGRA points. We find nearly 
perfect quantitative agreement between the expectation and MC data for 
wide parameter regions. In section 4, we show that systematics dues to 
the X2 velocity distribution could be corrected by the hardest lepton's Pt 
distribution. We also show expected errors on ino and slepton masses. For 
the most optimistic cases where the end point of the lepton invariant mass 
distribution of the three body decay m^^°'^^ is observed in addition to the edge 
of the mil distribution of the two body decay mf^°'^^, we can determine m^o , 
m^o and mj from those (almost) purely kinematical information. At least 
two degrees of freedom of the three mass parameters would be determined 
by our method if mj)^'" ^ 25 GeV. Section 5 is devoted to discussions. 

2 Distribution of lepton energy asymmetry 
with mil cut 

At hadron colliders, the second lightest neutrahno X2 would be produced in 
q and g decays, or in X1X2 P^i^ production. The decay X2 ~^ could be a 
dominant decay mode if it is open. Followed by / Ixi, the signal consists 
of a same flavor and opposite sign lepton pair associated with some missing 
momentum. It is one of the most promising SUSY signals at hadron colliders. 

The decay process X2 ~^ ^^^f ~^ Xi^i^^f is described by two body kine- 
matics and very simple. The mu distribution of the lepton pair from the X2 
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cascade decay is 

-— = (1) 

where 

" mj • 

The decay distribution is flat in mfi. The only physical information we 
can get from the mu distribution is therefore the value of the end point. It 
constrains one combination of the three masses involved in X2 decay, as one 
can see in Eq.(^. 

In the rest frame of the second lightest neutralino, the energy of h is a 
function of m^o and rrij, while Ei^ also depends on mu and fn^o: 

m|o - mf ml + mj - m?o 
Ei, = ^ ^, Ei, = — ^ (3) 

X.2 A2 

The angle On between the two leptons in the X2 ^^^^ frame is obtained by 
solving 

ml = 2Ei,Ei,{l-coseu) (4) 

^ = for mil = 0, while = vr for mu = m^^^. 

In Eq.(^), we see that Ei^ is monochromatic in the X2 ^^st frame. As 
a result, the energies of the two lepton are, asymmetric. The ratio of the 
transverse momenta of the leptons, which we call the transverse momentum 
asymmetry = Pt2/ Pti (-^ri > ^1:2)1 provides another information on the 
decay kinematics.^] However, Pip^ and Pip2 depend on the parent neutralino 
momentum, unlike the Lorentz invariant quantity mu. The X2 velocity dis- 
tribution in turn depends on and m^, although the X2 decay distribution 
in the X2 ^^st frame itself does not depend on them. 

-'^One may also use lepton energy ratio Eii/Ei2. In general, Pr distribution reflects 
sparticle masses much better than energy distribution. 
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This distribution has been used in global fits of MSUGRA parameters; 
Arp distribution "data" for one MSUGRA point generated by Monte Carlo 
simulator are compared to those of different MSUGRA points[^, ^. In this 
model, all sparticle masses depend on a few universal soft breaking parame- 
ters such as mo, M, tan /?, etc. When we compare different MSUGRA points, 
we therefore change both the parameters of the X2 decay, m^o, m^o and m^-, 
and the parameters of X2 momentum distributions rriq and rrig at the same 
time. Therefore it was considered to be less important compared to invariant 
mass distributions. 

However it is possible to make a more direct connection with the first set 
of mass parameters fn^o, m^o and rrij if a moderate mu cut is applied 0. 
When mil is small compared to 171"^^^, the lepton and anti-lepton nearly go 
in the same direction. Then the lepton momentum asymmetry becomes less 
sensitive to the parent neutralino velocity. Even after the smearing due to 
the boost of x?,, = Ei/E2\ma=o still can be extracted from the peak of 
At = Pt2/Pti Q; 

2 2 
m-o — mf 

(or IMD - A% = (5) 

mf — m-Q 

I Xi 

therefore A^"^^ constrains the mass parameters involved in X2 decay, just as 
^max ^Qgg^ Note A^ has monotonous dependence on all parameters while 
the mil edge might be accidentally insensitive on mj. 

In this paper, we study the power of the At distribution in the low mu 
region {mu < m^^^/2) to constrain the kinematics of the cascade decay 
X2 — ^ — ^ Xi^^- We chose the points shown in Table 1, but our method 
can be applied in generic MSSM studies. Unlike the common approach to 

^The lepton from X2 decay and I decay cannot be distinguished in the experiment. It 
is understood that yl^ means when A'^ exceeds one. 

■^We thank to M. Drees for pointing out this. 
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mo 


Ms 


m^o 


rriT 

l-R 


m^o 


^11 


A% 


IK[| 


100 


150 


135.5 


120.7 


65.2 


51.8 


0.368 


point 5|0, ^ 


100 


300 


233.0 


157.2 


121.5 


109.1 


0.336 


point 5-2 


115 


300 


233.2 


167.1 


121.6 


111.6 


0.496 


point 5-3 


120 


300 


233.3 


170.6 


121.6 


111.6 


0.565 


point 5-4 


125 


300 


233.3 


174.2 


121.6 


111.1 


0.646 



Table 1: Mass parameters and relevant sparticle masses in GeV for the points studied in 
this paper. ISAJET is used to generate sparticle masses. We also show corresponding 



mfl'^^ and in the table. 



immediately go into full MC simulations, we first study the decay distribution 
for fixed neutralino velocity (labeled by the boost factor 7^0 and the pseudo- 
rapidity r]^o) r{AT{j^o,rix'^)). 

The distribution we observe in experiments is expressed by convoluting 
the distribution with the velocity distribution of F{j,ri), as follows; 



da {A' 



T) 



(6) 



the measured distribution is also affected by cuts on I^t, Mes, etc. However 
it is still useful to know how r(ylj'(7^o, r^^o)) depends on the underlying mass 
parameters and the X2 velocity. 

In Fig. 1, we show the At distribution with/without invariant mass cuts 
and cuts. Here we take the IK point and = 1.9 , 77^0 = 0.2. The 
distribution is quite easily obtained by numerical integration. 

The distribution without upper mu cut has some structure around At = 
0.3 (top solid histogram), but it is insignificant. (Here we took the events 
with mil > 12 GeV because large backgrounds from virtual photons are 
expected for mu < 12 GeV [|n|). With the cut P| > 10 GeV and the same 

X2 velocity, events with At < 0.1 are hardly accepted, and the distribution 

"^The decay distribution F depends on 7 and ry through the Lorentz boost of all mo- 
menta, which is implicitly shown as ^7^(7, 77), or P^{'y,ri). 
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IK point 

0.05 
0.04 
0.03 
0.02 
0.01 
0.00 

0.2 0.4 0.6 0.8 1 

Figure 1: At distribution for the X2 momentum (7,?7)=(1.9, 1.2) without P| cut (soUd), 
and for PJj. > 10 GcV (dashed). The upper histograms are without upper mn cut while 
the lower histograms are distributions with mu < 25 GeV. 

is roughly flat between 0.2 < At < 0.3 (top dashed histogram). When the 
lepton energy in the X2 ^^^^ frame is small, the acceptance efficiency of the 
events strongly depends on the velocity of because of the Pj. cut. The 
At distribution would depend on the cuts and the distribution of X2 velocity 
introducing systematical errors to the analysis. 

On the other hand, once a moderate niu cut is applied, the decay distri- 
bution becomes nearly independent of Pj, cuts (bottom histograms). Here 
we integrate the region between 12 GeV< mu < 25 GeV ~ The 
distribution has a peak at At ~ A% — 0.368. The peak is outside the small 
At region affected by the Py cut. It is also clear from the plot that the shoul- 
der of the distribution without mu cut comes from the events with mu < 25 
GeV. Note that cos^^ = 0.86(0.44) for mu = 12(25) GeV in the X2 rest 
frame, therefore the angle between the lepton and the anti-lepton in the pair 
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is rather small with the uiu cut. 

To put it differently, events above m™^^/2 are merely backgrounds to the 
A'^ measurement. This is easily understood when we consider the lepton 
configuration near the mu end point. The two leptons go in exactly opposite 
directions and the asymmetry is modified maximally when one of the leptons 
goes to the direction of X2 momentum, A = / E}f° = A£;|mj;=mmax x 
where Ae = 0.29 and (3 = 0.855 for Fig. 1. The lepton energy asymmetry in 
the laboratory frame can range from nearly to 1 due to the boost. 

IK point 

0.80 



0.15 
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0.05 



0.00 

0.2 0.4 0.6 0.8 1 

At 



Figure 2: At distribution under different mu cuts. (7,??)= (1.4, 0.2), and P| > 10 GeV 
for solid and dotted histogram, while the dashed histogram is for 7 = 2.3 and t] = 0.2. 

It is worth noting that the At distribution peaks at smaller value of At 
as one increases the mu cut. In Fig. 2, we show distributions with different 
mu cuts, GeV< mu < 10 GeV (solid narrow), 10 GeV < mu < 20 GeV 
(dashed), 20 GeV< mu <28.3 GeV (solid wide), for 7^0 = 1.4 and r]^o = 0.2. 
The distribution has a sharp peak at a position consistent with A% for the 
sample with mu < 10 GeV. A^^^ = 0.323 for the same neutralino velocity for 



10GeV<m|,<20GeV 
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10 GeV < mil < 20 GeV (dashed histogram). This shift cannot be explained 
by Ae deviation from A%{Ae = 0.363(0.354)) for niu = 10(20)GeV, but it 
comes from the smearing of At distribution for finite lepton angle. 

The dotted histogram shows a distribution for a higher neutralino veloc- 
ity 7 = 2.3 and rj = 0.2 with 10 GeV < mu < 20 GeV. The peak position 
is shifted very little, A?,'^'' ~ 0.321, therefore it may still be used to deter- 
mine the decay kinematics^. On the other hand, the distribution off the peak 
depends more on the neutralino velocity. Using the whole distribution intro- 
duces a dependence on the X2 momentum distribution, and the fit would be 
more assumption dependent. 

The distribution is more and more smeared out and peaks at a lower At 
for larger mu cuts. The dependence on the X2 momentum is also bigger for 
the large mu sample; A^^""^ = 0.26 (0.24) for 7 = 1.4 (2.3) and 20 GeV < 
mu < 28.3 GeV. (Only the distribution for the former is shown in the fugure.) 
The distribution is shifted to smaller At reducing the acceptance of the 
Pt > 10 GeV cut. Some information on the neutralino velocity distribution 
is therefore necessary to deduce the neutralino decay kinematics from the At 
distribution while increasing the mu cut in order to increase the statistics and 
remove virtual photon backgrounds. This will be discussed in detail in section 
4. 

Note that m™^^ ~ 50GeV for IK, therefore requiring mu < 25 GeV re- 
duces the number of events in the sample by 1/4. The reward is a distribution 
which is less sensitive to Pt cuts and to the X2 velocity distribution, and a 
simple correspondence to the quantity in the X2 ^^st frame. 

Finally we demonstrate sensitivity of the At distribution to the slepton 

^Peaks are determined by fitting the distribution near the peak to a polynomial fitting 
function. 
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point 5 (ttiq— 1 OOGeV) 



point 5 (mo=115GeV) 




a) 




b) 



Fi gure 3: At distributions for different invariant mass cuts, j^^o =3.1, /3j^o = 0. The 
distribution for tight mu cuts are scaled by a factor of 5. 



mass. We first compare distributions with different slepton masses, P5 (mo = 
100 GeV) and P5-2 (mo = 115 GeV) in Fig. 3a) and 3b). Here we try a 
relatively large 7^0 in order to have a substantial effect from the X2 boost 
(7 = 3.1 1] = 0). Still, the distributions arc clearly peaked at At ~ 0.32 
(P5) 0.48 (P5-2) for events with mu < 50 GeV, while the distribution with 
mil > 50 GeV does not show any structure between At — 0.1 to 1. Note 
mil — 50 GeV roughly corresponds to half of m]^^^ again. 

In Fig. 4, we compare distributions with different mo. Peak positions 
shift from 0.3 to 0.63 as one changes mo by 25 GeV. If systematic errors 
are negligible and M is fixed, the sensitivity to mo would be Smo ~ 1.6 
GeV for SA — 0.02 (as will be found in section 3). The peaks are consistent 
with A% = 0.33 (for mo = lOOGeV), 0.49 (for mo = 115GeV) and 0.65 (for 
mo = 125 GeV). In section 3, we will find similar agreement for full MC 
simulation data, establishing the correspondence. 
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point 5 




Figure 4: At distributions for different slepton masses. Cuts are 12 GeV < mu < 50 
GeV, P| > 10 GeV. (7,7;) = (3.1,0) 

3 Monte Carlo simulations 



We are now ready to perform full Monte Carlo simulations to check the 
observations made in section 2. 

We use ISA JET 7.42||l0l to generate SUSY events. The generated events 



are analyzed by the simple detector simulator ATLFAST2.21 |]T2[. The cuts 
to remove the SM backgrounds down to a negligible level have already been 
studied in H; they are summarized as, 

IK (inclusive 3 lepton channel) §]; 

For this point, M^s and I^t cuts are not efficient because of the light g. A 
third tagging lepton from X2 or Xi decay is required. When three leptons are 
in a same flavor, the pair of lepton with smaller AR is selected as a lepton 
pair candidate. 

• Two opposite sign same flavor leptons with PJp > 15 GeV. 

• Third tagging lepton with Pj, > 15 GeV. 
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• Lepton isolation; No Pt > 2 GeV track within a AR < 0.3 cone 
centered on the lepton track. 

• ^T> 200 GeV. 
point 5 0; 

• 4 jets with Pti > 100 GeV and Pt2,3,4 > 50 GeV. 

• Meff = Pt,i + Pt,2 + Pt,3 + Pta + > 400 GeV. 

• > max(100 GeV,0.2Mefr). 

• Two isolated leptons with Pj, > 10 GeV, |?7| < 2.5. Isolation is defined 
as less than 10 GeV energy deposit within a AP < 0.2 cone centered 
on the lepton track. 

We generate 2 x 10^ events for each point. This roughly corresponds to 5 
for IK, and 100 fh^^ for point 5. We present distributions without cuts 
on Meff, jet Pt and Bfr- In previous simulations 0, 0; the acceptance 
is roughly constant for all value of mu, therefore those cuts are expected 
not to modify the lepton distributions substantially. Note that substantial 
acceptance for events with mu < m^^^ /2 is crucial for using the information 
from the At distribution, as we have seen in section 2. 
We keep lepton isolation cuts; 

• Less than 10 GeV (15 GeV for IK) energy deposit within a AP < 0.2 
cone centered on the lepton track. 

• No jet within a AP < 0.4 cone centered on a lepton track. 

^We use jet finding algorithm of ATLFAST. jet cone size is Ai?j < 0.4, Jet finding 
algoritlim requires 1.5 GeV of minimum energy deposit for the cluster seed, jet cone size 
Ai?j < 0.4, 10 GeV minimum total energy. A resulting cluster with energy more than 15 
GeV is called jet. 
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The acceptance of events turns out to be too high by factor of 3 for point 
5 compared to a full analysis including jet related cuts [§|, §, 01- This factor 
is taken into account when we interpret the fitting results.[] No plot or fit in 
this section contain SM background, while SUSY background is included. 
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Figure 5: (77,7) distribution of X2 for point 5. 



In the previous section, we have already seen that the At distribution is 
somewhat dependent on the parent neutralino velocity (7;^o, rj^o). In Fig. 5, 
we show the 7^0 and rj^o distribution for point IK. Here one can see that rj^o 
is roughly within \ri\ < 1. The X2 can be very relativistic; '-f^o could be much 
larger than 2. A modification of the At distribution due to Lorentz boosts 
is expected unless some mu cut is applied. The Pt distribution is shown 
in Fig. 6. Here we plot the distribution of higher (lower) of two lepton Pt 
for dotted(solid) line. The first (higher) lepton P^ can be a few times higher 
than its most probable value, refiecting the existence of relativistic X2 ^^e 
signal sample. 



^The number of the selected events for the IK point is 7000 between 10 GeV to 20 
GeV even for the small luminosity of 5 Therefore, the systematic errors would be 
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a) b) 
Figure 7: At distribution for a) IK and b) point 5. without upper mu cut. 
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We now study the asymmetry distribution in Fig. 7. The plot for Point 
IK (Fig. 7a) shows a smeared peak at At ~ 0.36, but the peak is rather fiat 
at the top. For point 5 (Fig. 7b), the distribution has even less structure, 
especially when mo > 115 GeV. Although global fits of the distributions 
must give us information on the neutralino decay kinematics, the power to 
constrain neutralino decay parameters would be limited if we try to analyze 
models without the constraint between soft breaking parameters. 



' J3 500 - 




a) b) 

Figure 8: At distribution with a 100 (GeV)^ < ml < 800 (GeV)^ cut 

In Fig. 8, we show Ax distributions with mu cut, 100 (GeV)^ < mfi < 
800 (GeV)'. We find a narrower peak for point IK compared to the case 
without invariant mass cut. For point 5, improvement of the signal distri- 
bution is clear. The peak position moves right as mo is increased, and it is 
consistent with Fig. 4. Note that for point 5 (IK), mfi < 800 (GeV)^ corre- 
sponds to cos6';; = 0.72 (0.29). The angle between the two leptons is smaller 
for point 5, which explains the substantial improvement for point 5. 

In Fig. 9, we show the distribution of the events with mu < 10 GeV 
for point IK. The peak is now nearly delta function like, and it agrees with 
dominant for this point. 
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Figure 9: At distribution for point IK. ma < 10 GeV 

A^. For point IK, the number of signal events in this range is statistically 
significant. If events in this mass range can be used, we should be able to 
make a direct A% measurement without any systematics. However, there 
could be a significant background for the events below mu < 12 GeV as 



recently discussed in [11 



We now fit the MC data to a phenomenological fitting function to de- 
termine the peak positions and the associated errors. The fitting function is 
chosen as follows: 



NiA) 
N{A) 



No exp -0.5 X 



A-Ao 



for A< Ac\ 



iVoexp(-/(A-Ao)) for A> 



(7) 



where parameters Aq, /, Nq and a are determined by minimizing using 
the program MINUIT0. 

Results of these fits are shown in Fig. 10. For IK (Fig. 10 a)) we fit 
the At distribution of events with 10 GeV < niu < 14.14 GeV and find 



^Here we take a completely phenomenological assumption for the fitting function, how- 
ever it is much better to use the fitting function based on the neutralino velocity distribu- 
tion calibrated by the first lepton distribution. See section 4. 
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c) d) 

Figure 10: Fits to the MC data using the fitting function 0. a) IK, b) point 5, c) point 
5-2 and d) point 5-4. 
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Aq = 0.3408 ±0.01. The peak position is smaller than = 0.36 (defined in 
Eq.(^). However, the distribution for fixed neutralino velocity {'~f,ri) = 
(1.4, 0.2) indeed peaks at 0.34 if 10 GeV < ma < 14.14 GeV, consistent with 
the full MC. As discussed earlier, the peak position does not depend strongly 
on the choice of the neutralino velocity when such tight ma cut is required; 
this explains the agreement. For point 5 (Fig. 10, b)-d)), Aq = 0.324 ±0.005, 
0.491 ±0.012, and 0.675±0.018 for mo = 100, 115 and 125 GeV, respectively. 
The mil cut dependence is rather small; A'^ is 0.33, 0.49, 0.65 respectively, 
already consistent with the fitj^. Recall that the total number of events is 
a factor 3 too large since we ignored jet related cuts. The corrected error 
for 100 fb-^ luminosity is 0.009, 0.02 and 0.03 for point 5, 5-2, and 5-4 
respectively, assuming statistical scaling. 

4 Model independent mass determination 

The second lightest neutralino might arise from squark and gluino decays 
at hadron colliders, therefore the X2 velocity distribution should depend on 
mq and mg. One may in principle fit the whole distribution to determine 
model parameters completely, but various systematic errors could prevent a 
complete understanding of the event structure. We wish to stay with the 
distribution which is less model dependent and free of systematics. Invariant 
mass distributions are a well established candidate for such a distribution. In 
the previous sections we argued that the peak position of the A^ distribution 
can be almost independent of the X2 velocity distribution if certain cuts are 
applied on mu. 

In this section, we will find that the remaining minor 7^0, rj^o distribution 

^Note that the peak positions are at larger At compared to Fig. 4. This is because we 
select the events below mu < 28.3 GeV for Fig. 8, while it is mu < 50 GcV in Fig. 4 . 
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dependence may be removed by looking into the first lepton Pp distribution. 
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first lepton 
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Figure 11: &) At and b) distributions for different X2 velocity (7,??) = 
(1.1, 0.2), (2.2,1) and (3.5,1.3). A^^''^ of each distribution is also indicated in the figure. 
We also show the histogram for my = 117.68 GeV in a) for comparison. 



In Fig. 11a), we show the At distributions for different X2 velocity. We 
took point IK and 12 GeV < mu < 25 GeV, therefore the distribution is 
somewhat dependent on the neutralino velocity especially when j^o is small. 
A^^^ shifts from 0.31 to 0.29 between the representative neutrahno velocity. 
In the same figure we also show distributions for mf = 117.68 GeV. In this 
case the peak moves from 0.42 to 0.37. The velocity dependence is slightly 
stronger than for the IK point. Although the peak position itself does not 
depend too much on the velocity, this certainly suggests some systematics 
would come into the fit to the decay parameters. 

The X2 velocity distribution strongly affects the hardest lepton distri- 
bution, as one can see in Fig. lib). Here the three distributions correspond- 
ing to Fig. 11a) have totally different end points. We can imagine that 
the systematics coming from the neutralino velocity dependence would be 
substantially reduced if the P^ distribution is included in the fit as well. 
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The At and distributions can be expressed as convolutions of the 
neutrahno velocity distribution and neutralino decay distributions as follows; 

a{AT) = y'rf7^orfr/^oF(7^o,r/^o) X r(Ar(7xO,^x§)) , (8) 
a{P^) = y'(i7^o(i7;^oF(7^o,r;^o) X r(P^(7^o,r;^o)) . (9) 

The neutralino velocity distribution F{'~f, rj) depends on parent sparticle masses, 
while the decay distributions in the laboratory frame T{At) and r(Pj,) de- 
pend on m^o, m^o and m^- implicitly. Various cuts would be applied to ex- 
perimental samples of events, therefore these equations are rather schematic. 
Note that the two distributions have different neutralino velocity dependence. 
The ?7^o and 'y^o dependence tend to cancel in AT{'y,r]), while the transverse 
momentum in the laboratory frame Py(7, 77) keeps increasing with 7. Hence 
a measurement of the Pj, distribution must be very useful for correcting the 
minor dependence of the At distribution on rj^o and ■j^o . 

The parent neutralino velocity can be decomposed into a boost 7^ from 
the X2 rest frame transverse to the beam direction, followed by a boost 71, 
in the beam direction. The At and Pj. distributions depend on the jt dis- 
tribution while the latter distribution has no effect on them. This can be 
seen in Fig. 12 a) and b). We show three At and Pf distributions, for 
X2 {ix-iVx) ~ (1-1) 0-2), (1-2,1) and (2.7,2). The 3 points have a common 
feature, 

Bt{i.v)^ ^ (10) 
£^/i(at X2 rest) 

The Pj, and At distributions of leptons are very similar as one can see in 
Fig. 12 a) and b). 

This observation is based on a numerical integration which now takes 
into account the cut \rii\ < 2.5, in addition to 12 GeV < mu < 25 GeV, and 
the Pj, > 10 GeV cut. The effect of the rii cut turns out to be very small. 
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a) b) 

Figure 12: At and distributions for different neutralino velocities but the same 
transverse boost factor Bt- Solid histogram is for {-j,r]) = (1.1,0.2) while bar graphs are 
for {-f,T]) ~ (1.2, 1) and (2.7,2) . Errors of numerical integrations are shown as bar size. 
Two bars in a same bin are almost coincide. A distribution with a different transverse 
boost factor is shown by the dotted histogram. 



We checked numerically that the distributions with common Pj. end points 
are roughly the same with these cuts. On the other hand, the At distri- 
bution has significant Bt dependence as one can see from the distributions 
for Bt — 1.7 (dotted histograms). This suggests that one only has to know 
the 7r distribution, which could be reconstructed from the Pj. distribution. 
Schematically, one can write 

aiAT) = JdBTF{BT)xr{ATiBTiv,l))), (11) 
a{P^) = ldBTF{BT)xr{PUBT{v,m- (12) 

T{At{Bt)), and T{Pt{{Bt)) are implicit functions of ino and slepton masses, 
and one can fit to experimental data to obtain those mass parameters in 
addition to F{Bt). Of course, one must also study the effect of ^t, Mes, 
and Ptj cuts and detailed MC simulations are necessary. 

Given the indication that the dependence on the X2 velocity distribution 
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can be corrected directly from the Pj. distribution, we now use the error on 
A^^^ and mu end points to determine fn^o, m^o, and mj. As we have seen 
in the previous sections, Aj^'^^ depends on the mu cut, but we assume the 
statistical uncertainty of A^^^ can be translated into that of A%; i.e., we 
assume that the correlation caused by only using events within a finite range 
of mil values would be small. 




20 40 60 80 100 



Figure 13: The mu distribution for point IK. 

We take the IK point as an example; the point is interesting because 
both the edge of the mu distribution due to the two body cascade decays 
^2body ^j^g point of the three body decay X2 ~^ ^^X?) can 
be seen. (See Fig. 13.) This is because the right handed slepton coupling to 
wino and higgsino is essentially zero, therefore the two body decay coupling 
is suppressed. The measurements of mf^°'^^ , Aj,, and m^^°'^^ = m^o — m^o 
are potentially sufficient to determine all sparticle masses involved in the X2 
cascade decay. 
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Figure 14: Contours of constant Ax^ = 1,4,9 for the IK point. We set Sm^a = 0, 
5m'f^°'^^ 0.5 GeV and (5A^ = 0.007. Only the contours near the input value are shown. 



In order to demonstrate the importance of the measurement of At, we 
first show the expected constraints on mj and m^o when m^o is fixed. We 
assume that and mfj^"'^^ are measured within errors of 0.007 and 0.5 GeV, 
respectively (Fig. 14). Ax^ is defined as 



(13) 



^^2body J 

Here, A^, 6A^^ and mu = mf^"'^^ and 6mfl^"'^^ are 'data' and the error, 
while, A% and m^^ are functions of the ino and slepton masses. The result- 
ing Ax^ = 1, 4, 9, . . . contours roughly correspond to la, 2a, 3(J, . . . errors on 
the parameters. The errors on m^- and m^o^ could be of the order of 1% or 
less, consistent with the previous fits in Note, however, that they did 
not identify the origin of the peak structure and used the whole At distri- 
bution for the fit. As we have stressed, this fit will depend on assumptions 
about parent squark and gluino masses, while our fit relies solely on the peak 
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position, directly constraining m^o, m^o and m^-. 



For the IK point, one can also determine m^^ ° ^ . The errors on the 



masses under the three constraints would be substantially larger than those 



shown in Fig. 14 (where m^o is fixed). This is due to correlations between 



the constraints. This can be seen in Fig. 15, where m^o and m^- are shown 



as functions of ° ^ for fixed values of and mfi ° ^ . Even in the limit 



where and mn ° ^ are known exactly, an error on mn ° ^ of the order of 



1 GeV would result 5 GeV errors on m^o and mr. 
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Figure 15: m^o and TOj as the function of m^^ ° ^ when A'^ = 0.37 and m^^ ° ^ = 51.9 
GeV. 

Assuming an error on m^ ^"'^^ , 6mfi = 1 GeV,^ m^o , m^o , and mj 
are constrained within ~ ±8 GeV, without assuming any relation between 
ino and slepton masses. The error is large compared to those expected from 

LC experiments, however it still makes an impressive case where sparticle 

"'^''Our assumptions of the errors for mu endpoints are substantially conservative to those 
found in hteraturc[^ Q 
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masses are determined without relying on model assumptions. |^ 

Note that the m^/mp, ratio would be constrained strongly. Assuming 
5At < 0.007, Sniec^i^i < 0.5 GeV, 6m^l'°'^^ = 4 GeV, we obtain 5(me/m^) = 
2.5 % for Ax^ < 1, and 7% for Ax^ < 9. 

Several comments are in order. The background in the region mu <^ 
j^max j^^g^ studied carefully. For example, SM ttU production could be 
important in the low mu region. Note that full amplitude level studies of Wj* 
production have been performed for the background process of X2 Xi^ 3/, 



and large background was found in the mu < 10 GeV region [|Tl|]. It has also 
been pointed out that T production is an important source of background 
when mu < 12 GeV. However it is unlikely that the background distribution 
has a peak at ^ 0. A peak of the signal distribution may be observed 
precisely on the top of such backgrounds, especially when signal rates are 
high enough to allow precision studies. Besides, one only needs to require 
mu < mfl^^/2 to see structure in the At distribution. The peak position 
that may deviate from A% could be corrected from the distribution in an 
almost model independent way. 

Recently, it was pointed out in [0] that one can obtain the same infor- 
mation by taking the ratio of the end points of the invariant masses of jet 
and lepton(s). Their analysis was carried out for point 5. The dominant 
cascade decay process is g — X2I followed by X2 ~^ / — ^ Xih- Jets 

from squark decays are substantially harder than the other jets, and can be 
identified. A correct set of a jet and a lepton pair originating from a squark 
decay is then selected by requiring that muj < 600 GeV for one of the two 

hardest jets j, and muj> > 600 GeV for the other jet j'. The end points of 

-'^^For point 5, end points of mu, miq, muq distributions in addition to the lower end 
point of muq distribution when mu > m^'^^ 12 is required determine to^o mass within 
0(10%) model independently. @] 
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the invariant mass distribution mi^g and mi^i^g are expressed as simple ana- 
lytical functions of rrig, mj, m^o^^^. One can reconstruct the rrii-^g end point 
by choosing the combination of the first lepton and the jetQ 

Although each end point is 10% to 4% smaller than expectation depending 
on jet definition, the ratio 



Hq 



^llq 



1 1 



(14) 



agrees with the expect at ion .P^ The fitted value ranges from 0.87 (Ai? = 0.4 
for jet definition) to 0.877 (Ai? = 0.7) while the expectation is 0.868. The 
range corresponds to A\ = 0.321 to 0.30 while the expectation is 0.327. Our 
fit gives Aq = 0.324±0.009 for the same point. The comparison of systematics 
might be an interesting topic for future studies. Our A^^^^ analysis may be 
performed even if jets and leptons in the same cascade decay chain can not 
be identified, therefore it can be applied in a wider context. 

It is also interesting to reconstruct the kinematics when both X2 ^ ^R^ 
and X2 ^ open and the branching ratios are of the same order. In 

addition to the two edges of the ma distribution m™'*^(low) and m™'^^(high), 
one should be able to observe two peaks in the A^ distribution Aj:^ and A^^ , 
corresponding to the two decay chains. By comparing At distributions for 
mil ~ rnii(\ow)/2 and mii(\ow)/2 < niu < mii(h.igh.)/2, one should be able 

^'^ Note that the efficient selection of the first lepton for mig distribution relies on large 
lepton energy asymmetry. However as ^ — s- 1, the end points of mi-^j and mi^j tend to 
coincide, therefor it may not be a problem. 



13 When m? - m|o > m|o - m?, mig/muq = y^(m? - m|o)/(m|o - m|o). 

Another potential problem of the analysis in Q is that rriq — m\^^^ — 145 GeV is 
almost as small as m^o = 122 GeV. Being the end point of the muj distribution requires 
the x\ from the decay chain to be very non-relativistic in the q rest frame. This should 
reduce I^t toward the end point. In general the muj and mij end points correspond to 
different kinematical configurations; attention must be paid to the consequence for relative 
efficiencies. 
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to determine proper sets of the mu edge and the peak, because the peak at 
Aj^^ can be hardly observed for rriu > mii{\ow)/2, while the peak at A^"* can 
still be seen. Note that there are four parameters for four constraints in this 
case, therefore one can in principle solve for all mass parameters. 

5 Discussion 

The second lightest neutralino X2 would be copiously produced from q and g 
decays at the LHC, and xtx2 production is an important mode for the Teva- 
tron. In this paper, we have studied the distribution of the Pj. asymmetry. 
At = Pt2/Pti^ of the lepton-anti-lepton pair that arises from the cascade 
decay X2 ^ ^ Xi^^- We have found that the At distribution shows a clear 
peak structure in a wide parameter region if mu < m™^/2 is required. The 
peak position is insensitive to the parent X2 velocity distribution, and in the 
limit of mil ~ 0, it is understood as A%^ the ratio of lepton and anti-lepton 
energy in the rest frame of j^. The ratio AP^ is a simple function of m^o^ mj 
and m^o. 

We have also performed MC simulations for several representative points. 
Values of the peak position obtained by fitting MC data agree with those 
for X2 with typical velocity. This follows from the insensitivity of the At 
distribution to the parent neutralino velocity. The typical velocity could be 
estimated easily by using the hardest lepton Pt distribution. Therefore the 
At peak can be used to constrain m^o , and m^o . By using the edge of the 
mil distribution in addition to the At distribution, one can determine two 
degrees of freedom of the three mass parameters involved in the X2 cascade 
decay. When the end point of mu distribution of the three body decay 
X2 ~^ Xi^^ can be measured simultaneously, one can determine the all mass 
parameters describing X2 cascade decays. The analysis is entirely based on 
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lepton distributions and does not rely on jet energy measurements. 

The reconstruction of the X2 momentum distribution is of some impor- 
tance for our analysis. The hardest lepton Pt distribution should allow us 
to study the X2 velocity distribution independently from the q, g mass de- 
termination. In fact, the measurement of this distribution may allow one to 
constrain the kinematics of squark and gluino production. 

The fit proposed in this paper is reasonably model independent compared 
to the previous fits using the entire At distribution without mu cuts. It is 
amazing to see that the distribution keeps the information on the cascade 
decay kinematics. (Compare Fig. 7 and 10). The analysis can be extended 
to all cascade decays involving leptons, such as the gauge mediated scenario 
with NLSP slepton.[^, ^ The determination of the Aj- peak position is not 
disturbed even in the case where several sleptons contribute to signal lepton 
pairs. 

Note that model independent constraints on weakly interacting sparticle 
masses may be used to directly constrain the relic mass density of LSPs in 
our Universe. The density of such Big Bang relics is roughly proportional to 
the inverse of the pair annihilation cross section of the lightest neutralino. In 



the MSUGRA model, 1/cr ~ mj/m-o in the bino dominant limit. Il4| If the 
overall sparticle scale is constrained within 10%, an upper bound on the mass 
density could be derived within 20%. The improved determination of SUSY 
parameters at the LHC combined with improved astronomical observations 
might significantly constrain the remaining MSSM parameters. 

In this paper, we did not perform any MC simulation for Tevatron exper- 
iments. There the cleanest discovery process is the 3 leptons and missing I^t 
channel of X1X2 production and decay. It is possible to perform a parallel 
analysis to the one presented in this paper. However if mf^^^ is small (which 
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is likely due to the lower bound on mj of nearly 100 GeV), the number of 
events that satisfy 12 GeV < niu < m™'^^/2 would be small, where the lower 
ma cut is needed to avoid 7* and T backgrounds. 

The branching ratio of the mode Xo, ~^ could be small if other modes 
such us X2 ~^ Z,h... dominate. The decay X2 ~^ ^'^ ^e only two body 
decay channel accessible in MSUGRA model due to f mixing. The analysis 
would be substantially more difficult for this CclSG, clS T decays further into a 
jet or a lepton.[0] Selecting two tau leptons which go roughly into the same 
direction (small AR) should effectively work as an cut in our analysis. 
However, the At distribution of tau jet would be substantially smeared by 
the tau decay. 

When all two body decay modes are closed, the decay X2 Xi^^ often 
has a sizable branching ratio. The precision study of the three body decay 
distribution has been discussed in [T^. The ma distribution and the At 



distribution in the small mu region would give us information on neutralino 
mixing and on mr 

It would be interesting to check if our analysis can be extended to other 
cascade decays involving photons or jets Note that in the gauge mediated 
model with x1 NLSP, the decay chain X2 ~^ Xi^^ be associated with a 
photon from Xi ~^ G'y [|l^. Cascade decays involving a jet and a lepton 
or two jets may also be used for an asymmetry analysis, but selecting the 
proper combination of jets would be challenging. 
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